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SUMMARY 

I he reliability of scanning laser acoustic microscopy t,SLAM) for detecting 
surface voids In structural ceramic test specimens was statistically evaluated. 
Specimens of sintered silicon nitride and sintered silicon carbide, seeded with 
surface voids, were examined by SLAM at an ultrasonic frequency of 100 MHz In 
the as-fired condition and after surface polishing. It was observed that pol- 
ishing substantially Increased void detectability. Voids as small as 100 vim In 
diameter were detected In polished specimens with 0.90 probability at a 0.95 
confidence level. In addition, Inspection times were reduced up to a factor of 
10 after polishing. The applicability of the SLAM technique for detection of 
naturally occurring flaws of similar dimensions to the seeded voids Is dis- 
cussed. A Fortran program listing Is given for calculating and plotting flaw 
detection statistics. 


INTRODUCTION 

Silicon nitride (S1 3 N 4 ) and silicon carbide (SIC) ceramics are under 
Investigation as candidate materials for hot-section components In advanced 
heat engines (refs. 1 to 5). Because these ceramics can withstand higher oper- 
ating temperatures than their metallic counterparts, their use would result In 
significantly Increased fuel efficiency. Presently, state-of-the-art struc- 
tural ceramics exhibit wide variability In strength and low fracture toughness. 
These undesirable properties are generally attributed to flaws Introduced dur- 
ing fabrication processes In the form of voids, microcracks, and foreign mate- 
rial Inclusions (refs. 6 to 12). Flaws as small as 10 ym have been defined as 
critical; that Is, potentially failure causing (refs. 6, 8, and 9). 

Sensitive, reliable, nondestructive evaluation (NDE) techniques are needed 
to detect flaws In structural ceramics and reject parts containing critical 
flaws or concentrated flaw populations (refs 7, 8, and 13). NDE techniques can 
also aid In process optimization by Identifying the stages of fabrication dur- 
ing which flaws are Introduced (refs. 7, 13, and 14). Scanning laser acoustic 
microscopy (SLAM) Is an attractive NDE technique because of Its ability to 
Image surface and subsurface microflaws In real time. It Is applicable to 
denslfled ceramics, and has previously been shown to be capable of detecting 
. critical flaws In S1 3 N 4 and SIC specimens (refs. 5, 7, 9, 10, 15, and 16). 

To date, however, a complete statistically based evaluation of the reliability 


of SLAM for detecting failure-causing flaws In structural ceramics has not been 
accomplished. 

This report describes a study that was conducted to evaluate the relia- 
bility of SLAM for detecting surface voids In sintered SI 3 N 4 and sintered SIC 
specimens. The approach was to determine detection reliability for statisti- 
cally significant populations of seeded surface voids In specially prepared 
laboratory specimens. The effects of specimen thickness and surface condition 
on void detectability were Investigated. The applicability of the reliability 
results obtained for the seeded surface voids to naturally occurring Internal 
and surface-connected flaws Is discussed. 


STATISTICAL RELIABILITY THEORY 

The reliability of an NOE Inspection technique Is a quantitative measure 
of the ability of that technique to detect flaws of a specific type and size 
In a particular material. Reliability assessment Is probabilistic In nr.ture 
because Inspection results are Influenced by many variables. These variables 
Include flaw shape and orientation, material surface texture and microstructure, 
and equipment and operator performance (ref. 17). Methods for analyzing the 
reliability of NDE Inspection techniques are discussed In reference 18. This 
study was based on specimens containing seeded surface voids where the total 
number of seeded voids and their locations were known. Since an existing void 
was either detected or not detected (only two outcomes possible), SLAM relia- 
bility was determined by using binomial distribution statistics (refs. 17 and 
19). 


By using binomial distribution statistics, an Initial estimate for the 
true (unknown) probability of detection of voids of diameter d can be taken as 
(ref. 17) 


p-jj (i) 

where p - Is defined as a point estimate of true probability, S Is the number 
of detected seeded surface voids of diameter d, and N Is the total number of_ 
seeded surface voids of diameter d. There Is an uncertainty associated with p" 
because It Is calculated vor a relatively small number of Inspections. There- 
fore, a conservative confidence level estimate of the true probability Is pre- 
ferred. This estimate Is defined as the lower-bound probability p^. The 
lower-bound probability Is considered an appropriate measure of the reliability 
of an NDE Inspection technique (ref. 17) and Is used In this study to describe 
the reliability of SLAM. A lower-bound probability (of detection) p^ can be 
calculated from the following expression (ref. 17): 

,rx), R> x o-p/-* <2) 

where G is the selected confidence level. 

A statistically significant probability of detection Is 0.90 at a 0,95 
confidence level (ref. 18). It Is not sufficient just to have a high ratio of 
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voids detected to voids seeded to obtain 0.90 probability of detection at a 
0.95 confidence level using equation (2). Probability Is also dependent 
on the quantity of voids seeded. For example, If 10 voids 100 urn In diameter 
are seeded and all are detected, the probability of detection of 100 urn voids 
Is only 0.74 at a confidence level of 0.95. It Is necessary to have 29 voids 
detected out of 29 voids seeded to obtain a probability of detection of 0.90. 
Further, 0.90 probability of detection at a 0.95 confidence level means that 
there Is a 0.05 (1 - 6 ) probability that 0.90 Is an overestimate of the true 
probability of detection 


MATERIALS AND PROCEDURES 

In this study, specially prepared ceramic specimens seeded with surface 
voids were used to characterize the reliability of 100 MHz SLAM for detecting 
flaws In structural ceramic test specimens. Seeded surface voids (as opposed 
to another flaw type) were used because they could be easily Identified (number 
and location) and accurately characterized (size and shape) by visual methods. 
This allowed the Investigator to gather accurate flaw detection data during 
SLAM Inspections. The specimens and voids were characterized using surface 
profiling Instrumentation, metallography techniques, and optical and electron 
microscopy. SLAM Inspection of the specimens In the as-fired condition and 
after polishing was performed, and data on void detectability were gathered. 

The Inspection data were grouped according to void diameter and analyzed. 

Curves of probability of detection as a function of void diameter were 
generated . 


Ceramic Specimen Development 

Test specimens similar In composition to typical sintered SI 3 N 4 and sin- 
tered SIC were fabricated with seeded surfau. ids. The starting materials 
were -100-mesh SI 3 N 4 powder containing Y 2 O 3 and SIOj sintering aids and -100- 
mesh alpha-SIC powder containing boron and carbonaceous resin binders. A 
selected amount of either powder was pressed In a double-action tungsten car- 
bide (WC)-llned die at approximately 60 MPa to form a modulus-of-rupture (M0H) 
test bar. While still In the die, the bar was carefully dusted using a 
moisture-free aeroduster In order to remove excess powder from the top surface. 
Approximately 20 styrene dlvlnylbenzene microspheres of the same size (115, 80, 
or 50 pm In diameter) were positioned on the top surface of the bar along the 
longitudinal axis. The test bar was then pressed at approximately 120 MPa, 
removed from the die, and once again dusted. With the microspheres now 
Impressed Into the test specimen surface, the specimen was vacuum sealed In 
thin-walled latex tubing and cold Isopressed at 420 MPa. 

The procedure was repeated to form all the green test specimens. Differ- 
ent bar thicknesses were obtained by using different amounts of powder. The 
seeded specimens were then heated to 525 °C under a vacuum and baked for 45 min 
to allow the polymer microspheres to decompose. A crater, or surface void, was 
left at each position on the specimen surface where a microsphere had been 
Impressed. Following the vacuum heat treatment, the surface of each specimen 
was dusted to remove any remaining debris from the Intentionally created sur- 
face voids (craters). Finally, the SI 3 N 4 test specimens were sintered at 
2140 °C under 5 MPa of nitrogen pressure for 2 hr. The SIC test specimens were 
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sintered at 2200 °C under 0.1 MPa of argon pressure for 0.5 hr. The as-fired 
specimens were approximately 30 mm in length and 6 mm In width, while the 
thicknesses varied from approximately 2 to 4 mm. 

After Initial SLAM Inspection of the specimens In the as-fired condition, 
the specimens were surface treated In the following manner: 23 sintered sil- 

icon nitride (SSN) specimens and 2 sintered silicon carbide (SSC) specimens 
were Individually hand polished. The seeded void surface of each specimen was 
pressed against either 600-grit (for SSN) or 320-grit (for SSC) silicon carbide 
grinding paper attached to a rotating metallographlc polishing wheel. The 
opposite surface of the specimen was not specially prepared In any way. The 
surface of a single SSN specimen was diamond ground using 400-grit diamond 
grinding paper. After surface treatment, all specimens were reinspected using 
SLAM. 


Specimen and Void Characterization Techniques 

The surface condition of several representative specimens of each material 
In the as-fired condition and after surface polishing was evaluated using a 
surface profile measuring system. The surface condition of the diamond-ground 
SSN specimen was also measured perpendicular to the grinding direction. The 
surface profiler used a diamond stylus, 12.5 urn In diameter, In contact with 
the surface of the specimen to measure the peak-to-valley roughness. A 2 mm 
length of the surface of each specimen was profiled. The maximum peak-to- 
valley roughness was obtained from each profile. To obtain Information about 
the material microstructure, polished, unetched cross sections of representa- 
tive SSN and SSC specimens were examined at a magnification of 200 by using a 
metal lograph . Optical micrographs were acquired at different cross-sectional 
locations . 

A photograph of the seeded void surface of each specimen was obtained at 
a magnification of 5 to aid In locating the voids during optical, electron, 
and acoustic microscopy. The location of each void was determined co the near- 
est 0.2 mm relative to an x-/ coordinate system. Each void was then examined 
at a magnification of 200 b' using the metallograph. The void diameter was 
measured using a calibrates reticle Inserted Into the metallograph. The depth 
of the void was measured by focusing first on the surface of the specimen near 
the void and then on the bottom of the void. Subtraction of one reading (of 
the calibrated graduations on the fine focus control) from the other gave the 
depth of the void. Optical micrographs of representative voids for each spec- 
imen were obtained. In order to Investigate void morphology In. more detail, 
several specimens were examined at a magnification of 400 to 450 using a scan- 
ning electron microscope (SEM). SEM micrographs of representative voids were 
acquired. Voids were characterized In as-fired specimens and after surface 
polishing. 


Scanning Laser Acoustic Microscopy 

Figure 1(a) shows a schematic diagram of the operation of the scanning 
laser acoustic microscope. The SLAM makes use of a laser to detect mechanical 
distortions (of the order of angstroms) produced on the surface of a specimen 
by piezotransducer-generated, high-frequency ultrasonic waves traveling through 
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the specimen (ref. 9). In this manner, an acoustic picture of the specimen, 
Including surface and subsurface defect., such as voids, Inclusions, and cracks, 
Is obtained end displayed on a video monitor. 

Figure 1(b) Illustrates the experimental setup used when Inspecting 
ceramic specimens on the SIAM. The specimen Is placed on the SLAM stage over 
the piezoelectric transducer. A clear plastic coversllp, coated on one side 
with a thin (approximately 0.1 ym thick) film of gold, Is placed on top of the 
specimen. The purpose of the metallized coversllp Is to provide a mlrrorllke 
reflective surface for the laser. Typically, the surface of the ceramic spec- 
imen Is too rough to reflect light In a mlrrorllke manner. Distilled water Is 
used as a couplant between transducer and specimen and between specimen and 
coversllp. The transducer, located In a small well 0.5 mm below the stage sur- 
face, radiates continuous 100 MMz ultrasonic waves toward the specimen at an 
Incident angle of 10°. The longitudinal ultrasonic waves are transmitted 
through the water couplant to the specimen surface, where they are refracted 
(based on Snell's law). Only the shear wave component, traveling at an angle 
of approximately 45° from the vertical, Is utilized fo- SLAM Inspection of SSN 
and SSC. The Interaction of the shear waves with the top surface of the spec- 
imen sets up a ripple pattern on the top surface. The ripple Is transmitted 
through the water couplant to the coversllp, where It sets up a corresponding 
ripple pattern on the gold film. The peaks of the ripple vary In amplitude 
according to the Intensity of the ultrasonic waves producing them. 

A laser beam constantly raster scans an approximately 2.0 by 2.3 mm area 
of the coversllp, Is angularly modulated by the peaks of the ripple pattern, 
and Is reflected to a photodetector. The ; notodetector converts the modulated 
laser light to an electronic signal. This signal Is processed and used to 
create a real-time, black-and-white acoustic Image that Is displayed on a video 
monitor at a magnification of approximately 100. Generally, the brightest 
regions on the acoustic Image represent areas of high acoustic transmission 
through the ceramic specimen while the darker regions correspond to areas of 
low or no acoustic transmission. 

SLAM Inspection of a specimen was performed with the surface containing 
the seeded voids nearest the laser. The specimen was positioned such that the 
specific void to be detected was located near the center of the laser spot. 

This procedure reduced the possibility of acoustic Images of seeded voids being 
confused with acoustic Images of naturally occurring flaws which were similar 
In appearance. The x-y coordinates of the voids, obtained from optical photo- 
graphs of the specimen, were utilized In this procedure. Some specimens had to 
be rotated between 0° and 180® about the laser axis to obtain the best acoustic 
Image. Some specimens were slightly warped, and the tilt controls on the SLAM 
stage had to be constantly adjusted to maintain the optimum acoustic Image. 
Detection of a void was defined as the ability to discriminate the void from 
background noise (ref. 14) due to naturally occurring flaws and material arti- 
facts. The time It took to Inspect a specimen was noted and used as an Index 
of difficulty In finding the seeded surface voids In that particular specimen. 
Acoustic micrographs of representative voids In each specimen were obtained. 
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Data Grouping and Analysis 


At most discrete void diameters, an Inadequate number of voids were exper- 
imentally produced to provide a valid statistical sample. Also, there was an 
error Involved In measuring the actual void diameter. It would have been mis- 
leading, therefore, to calculate probability of detection for discrete void 
diameters. Instead, the voids were grouped Into small (10 pm) Intervals 
according to diameter, and the probability of detection was calculated over 
these Intervals. The optlmlzed-probablllty method (ref. 17) was used to fur- 
ther arrange the voids, because the number of voids In many 10 -pm Intervals 
wr,s still Insufficient for a valid statistical sample. This method Increases 
the size of the sample used to calculate probability by Including Inspection 
data from Intervals containing smaller flaws. I + s use Is Justified by assuming 
that the probability of detection Increases with Increasing void diameter 
(ref. 17). Appendix A Illustrates the use of the optlmlzed-probablllty method 
to obtain the probability of detection results In this study. Appendix B lists 
a Fortran computer program for grouping the void detectability data, performing 
probability calculations, and plotting the results In the form t* probability 
of detection (at a 0.95 confidence level) as a function of void diameter for 
each experimental data set. 


RESULTS ANO DISCUSSION 
Specimen and Void Characterization 

Table I Indicates specimen surface condition, density, and size for each 
material. The specimens were similar In surface condition and density to typ- 
ical sintered SI 3 N 4 (SSN) and sintered SIC (SSC) MOR bars. The specimens 
were grouped Into three discrete thickness ranges for SSN and two discrete 
thickness ranges for SSC. Throughout the remainder of this report, the SSN 
specimen thickness ranges 2.1 to 2.2, 2.9 to 3.2, and 3.8 to 4.1 mm are 
referred to as 2, 3, and 4 mm thicknesses, respectively. The SSC specimen 
thickness ranges 2.8 to 3.0 and 3.7 to 4.0 mm are referred to as 3 and 4 mm 
thicknesses, respectively. 

Figures 2 (a) to (e) show center and edge optical micrographs of unetched 
cross sections Illustrating the variation of porosity within typical specimens. 
The pores are the dark spots In the micrographs. The porosity was uniform for 
SSC but varied from center to edge for SSN specimens. Porosity was greater at 
the center than near the edge for SSN. 

Figure 3 shows an optical photograph of the seeded void surface of a 
typical test specimen. All the voids could be readily located from such photo- 
graphs. Figure 4 shows optical micrographs of a typical seeded surface void 
In a test specimen In the as-fired condition and after polishing. In as-fired 
specimens, the error in the measurement of the void diameter using the metallo- 
graph was estimated at only ±5 pm because the void perimeter was essentially 
circular and readily definable. A1 ' sr polishing, the thickness of each bar 
was reduced by approximately 10 to 30 pm and the measured void diameter was 
Increased by approximately 10 lo 30 pm because of spalling at the void perim- 
eter. The error In the diameter measurement after polishing was estimated at 
±10 pm. In this case, void borders were Jagged and hence more difficult to 
define than those In the as-fired specimens. 
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The diameter ranges of the seeded voids for each material, specimen thick- 
ness, and surface condition are shown In table I. Overall, the voids ranged 
from 40 to 165 yrn In diameter. Figures 5(a) to (c) Illustrate the typical 
morphology of the seeded surface voids. Figure 5(a) shows an electron micro- 
graph of a typical seeded void In an as-fired specimen. Figures 5(b) and (c) 
are side-view diagrams Illustrating the difference In shape between voids In 
SSN specimens and those In SSC specimens. From the diameter and depth data and 
electron micrographs,, It was Inferred that the voids were ellipsoidal rather 
than spherical. The voids In SSC specimens were generally shallower than those 
In SSN specimens. 


SIAM Inspection of Specimens 

Figures 6 and 7 show representative surface profiles and associated acous- 
tic micrographs for SSN and SSC specimens, respectively. In. the as-fired con- 
dition and after surface treatment. The acoustic micrographs (Including the 
appearance of the detected voids) are typical for all specimen thicknesses 
examined. Ultrasonic wave Interaction with the relatively large and random 
surface roughness typical of the as-fired specimen created small dark rings of 
Interference that mottled the acoustic Image for SSN or SSC specimens In the 
as-fired condition (see figs. 6(a) and (b) and 7(a) and < b) ) . The mottling 
masked the seeded surface voids making them difficult, If not Impossible, to 
detect. Regardless of specimen thickness, the acoustic Images of SSN and SSC 
specimens Improved considerably after polishing only the surface containing 
the seeded voids (compare figs. 6(b) and (d) and 7(b) and (d)). Void detect- 
ability was substantially enhanced since only negligible Interference was 
generated by the relatively small surface roughness typical of the polished 
specimens (see figs. 6(c) and 7(c)). (The residual surface roughness of the 
polished bars was attributed to numerous tiny Imperfections on the surface as a 
result of Incomplete polishing. These Imperfections are evident In the acous- 
tic Image of the polished SSN specimen (see fig. 6(d)) and are an indication 
of the sensitivity of 100 MHz SLAM. The optimally polished areas of the speci- 
mens had a mlrrorllke surface finish.) Many of the voids not detected during 
Inspection of the as-fired specimens were detected after polishing, while those 
originally detected In the as-fired specimens were much more easily detected 
after polishing. Individual specimen Inspection times were up to an order of 
magnitude shorter after polishing because of the Improved void detectability. 
Typical Inspection times for as-fired and polished specimens were 50 min and 
5 min, respectively. 

Figure 6(e) shows the surface profile of the diamond-ground SSN specimen 
taken perpendicular to the grinding direction. The diamond-ground surface Is 
visibly smoother than the as-fired surface but contains periodic roughness 
greater than the random roughness In the hand-polished surfaces (compare 
figs. 6(c) and (e)). The acoustic Image of the diamond-ground specimen, how- 
ever, Is similar In clarity to that of the hand-polished specimen despite the 
presence of strlatlons produced by grinding marks (compare figs. 6(d) and (f)). 
It Is believed that the latter result Is due to the fact that the roughness 
(grinding marks) of the diamond-ground specimen Is ordered and uniform as well 
as low (as compared to that of the as-fired surface). It Is obvious from the 
previously mentioned figures that void detectability Is directly related to 
acoustic Image quality. Therefore, results similar to those obtained for the 
hand-polished bars might be expected for diamond-ground MOR bars which are used 
In structural ceramic development and evaluation. 
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It should be noted that, In general, the acoustic Image of as-fired SSC 
specimens were slightly clearer than that of as-fired SSN specimens (compare 
figs. 6(b) and 7(b)). Hence, the voids were generally easier to detect In as- 
fired SSC than In as-fired SSN. The reason for this Is unclear, but It Is 
likely to be due to differences In surface topology and specimen microstructure 
between the as-fired materials. Although distinct topological differences were 
not apparent from the profile measurements, It Is possible that microscopic 
differences could exist on the surface of the materials. Moreover, It Is felt 
that the SSN, because of Its variable pore distribution (see figs. 2(b) and 
(v.)), may have been a more attenuating material than the SSC. A difference In 
acoustic Image with specimen thickness was noted only for the as-fired SSN 
specimens. In this case, the acoustic Images of the 2 and 3 mm thick specimens 
appeared slightly clearer than the Images of the 4 mm thick specimens. 


Reliability of Void Detection 

The seeded void detectability data obta'r.od using SIAM Is presented In 
table I. These data are presented In raw form for each material, specimen 
thickness, and surface condition in the void size distribution plots of 
appendix C. Figures 8 to 10, derived from the raw data given In appendix C, 
show the statistical reliability of SLAM In the form of probability of detec- 
tion (at a 0.95 confidence level) as a function of void diameter (POD) curves. 

It Is Important to note two competing factors that may have biased the POD 
curves. First, the locations of all seeded voids were known, and each test 
specimen was examined carefully at these locations. Though essential to estab- 
lishing probability of detection statistics, this procedure biased the detec- 
tion probabilities toward hlgher-than-normal values since neither the number 
nor the location of flaws Is known during normal part Inspections. A second 
factor, on the other hand, Is the fact that although the optlmlzed-probablllty 
method of data grouping Increases the size of the sample used In calcualtlng a 
probability value, the value obtained Is conservative because It Includes data 
from smaller voids. Each point on the POD curves corresponds to the largest 
void diameter contained In the Interval over which probability was calculated. 
Although a POD curve can be plotted through the midpoints of the void diameter 
Intervals, the conservative approach of reference 17 Is preferred to minimize 
the bias resulting from prior knowledge of void location. 

Figures 8(a) to (c) show POD curves for the seeded voids In as-fired and 
polished SSN specimens of thickness 4, 3, and 2 mm, respectively. These 
figures Illustrate the Improvement In void detectability after polishing. For 
the as-fired specimens of any thickness, no voids of any diameter were detected 
with 0.90 probability. Overall, only 45 percent of the voids examined In the 
as-fired SSN specimens were detected (129 voids detected out of 284 voids 
examined). After polishing, voids as small as 100 to 150 \im In diameter were 
detected with at least 0.90 probability. (The smallest void diameter at which 
0.90 probability of detection occurred was different for each specimen thick- 
ness.) Had a sufficient number of voids (large enough statistical sample) been 
seeded In the 50 to 100 »im diameter range, it Is felt that voids In this range 
would have been detected with at least 0.90 probability as well. Overall, 

97 percent of the surface voids examined In the polished SSN specimens were 
detected (262 voids detected out of 270 voids examined). Some of the SSN test 
specimens were slightly warped making uniform polishing difficult. As a 
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result, small areas or i. ( e surface of a warped bar remained with the as-fired 
roughness. The acoustic Image In these areas was mottled and void detection 
was difficult. Voids undetected as a result of the mottling caused the curve 
labeled “polished" In figure 8(a) to dip below 0.90 probability at void diam- 
eters between 85 and 125 pm. If the polishing had been uniform, It Is likely 
that all of the seeded voids would have been detected. 

Figures 9(a) and (b) show Pfcf) curves for the seeded voids In as-fired SSC 
specimens for the 4 and 2 mm spar Imen thicknesses, respectively. Voids as 
small as 100 to 150 pin In diameter were detected with a least 0.90 probability. 
(The smallest void diameter at which 0.90 probability of detection occurred 
was different for each specimen thickness.) Overall, 88 percent of the voids 
examined In the as-fired SSC specimens were detected (306 voids detected out 
of 346 voids examined). 

As discussed previously, the acoustic Image of as-fired SSC specimens was 
slightly less mottled than that of as-fired SSN specimens. Hence, better POD 
results were obtained for as-fired SSC than for as-fired SSN. Detection 
voids was still difficult, however, and Inspections for Individual as-fired SSC 
specimens were approximately as lengthy as those for as-fired SSN specimens 
(up to 60 min). If the Inspection time for Individual as-fired SSC specimens 
had been limited to that recorded for Individual polished specimens, the POO 
results would have been much poorer. Only representative SSC specimens were 
polished because polishing would not significantly Improve POD results for SSC 
specimens (most voids In the as-fired SSC specimens were found). Therefore, 
statistical data are not available for polished SSC specimens. To Illustrate 
typical Inspection results for polished SSC, a SSC specimen seeded with 13 sur- 
face voids was Inspected In the as-fired condition and after polishing. In the 
as-fired condition, 10 voids were detected out of 13 examined In an Inspection 
taking 45 min. After polishing, 13 voids were detected out of 13 examined In 
an Inspection taking 5 min. 

Figures 10(a) to (c) show the POO curves of figures 8 and 9 organized to 
Illustrate the effect of specimen thickness on void detectability for each 
material. Since a difference In the acoustic Imago with specimen th 4 ckness was 
noted only for the as-fired SSN specimens (clearer acoustic Image for thinner 
specimens), POD noticeably varied with thickness only In this case (compare 
figs. 10(a) to (c)). As shown In figure 10(a), better POD results were 
obtained the thinner the as-fired SSN specimen. POD did not appear to be a 
function of specimen thickness for polished SSN or as-fired SSC specimens (for 
the thicknesses Investigated) since the respective acoustic Image of each mate- 
rial did not change noticeably with thickness. If similar void diameter dis- 
tributions had been available for each specimen thickness for the polished SSN 
and as-fired SSC, the POD curves within figures 10(b) and (c) would be expected 
to be nearly Identical. 

The detectability results presented for seeded voids In SSN and SSC spec- 
imens Indicate that reliability of void detection and ease of Inspection are 
directly related to acoustic Image quality. Since acoustic Image quality Is 
highly dependent on specimen surface condition, present as-fired specimens of 
SSN and SSC need an Improved surface condition If they are to be Inspected for 
flaws In a time-efficient and reliable manner using SLAM. 
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In this study, only artificially seeded voids In polished SSN specimens 
were detected with 0.90 probability at a 0.95 confidence level. It Is 
believed, however, that this result can be extended to naturally occurring sur- 
face and near-surface flaws In smoothly ground structural ceramic test speci- 
mens. (For purposes of discussion, a near-surface flaw Is arbitrarily defined 
as one that Is not surface connected but Is wholly contained within two lengths 
of the surface.) To Illustrate, figure 11 shows acoustic Images of three nat- 
urally occurring flaws In a diamond-ground silicon nitride specimen (ref. 7). 
Flaw A Is a surface-connected Inclusion, and flaw B Is a surface-connected 
crack-like defect. Flaw C Is a near-surface defect having a crack-like appear- 
ance. Each flaw Is near 100 pm In either length or width. As was the case 
for almost all seeded voids In polished SSN and SSC specimens, flaws A, 3, and 
C were readily detected because of the clarity of the acoustic Image. Hence, 
flaws In this size range In diamond-ground structural ceramic test specimens 
would be expected to be detected with C' 90 probability at a 0.95 confidence 
level. It Is worth .toting that while SLAM Imaged all three flaws In figure 
11, flaw B could not be resolved by radiographic methods and flaw C could not 
be resolved by radiographic or optical methods. 


CONCLUSIONS 

Scanning laser acoustic microscopy was determined to be a time-efficient, 
statistically reliable technique for detecting surface-connected voids In 
structural ceramic specimens with smooth surfaces. Surface voids as small as 
100 urn In diameter In polished sintered silicon nitride specimens were detected 
with 0.90 probability at a 0.95 confidence level. Similar detection reliabil- 
ities were not achieved for voids in as-fired sintered silicon nitride speci- 
mens that exhibit rough surfaces. Additionally, Inspection time was reduced 
up to a factor of 10 after as-fired surfaces were polished. Seeded surface 
voids were used to establish reliability of detection statistics In this Inves- 
tigation because they could be accurately characterized In terms of their size 
and shape using visual techniques. Evidence was presented, however, showing 
that the detectability of surface-connected and near-surface flaws In specimens 
with diamond-ground surfaces may be similar to the detectability obtained for 
the voids In the polished specimens. Hence, It was Inferred that the detection 
reliabilities reported herein are applicable to surface and near-surface flaws 
In smoothly ground specimens. 
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APPENDIX A 


9PTIMI ZED-PROBABILITY METHOD OF DATA GROUPING 

Table II shows a sample set of Inspection data arranged In (approximately) 
10 ym intervals. The Intervals were set so that each void appeared In only one 
Interval. The steps of the optlmlzed-probabl 1 1 ty method of data grouping to 
obtain probability of detection results from this data are Illustrated In 
figure 12(a) . The results are displayed as a plot of probability of detection 
(at a 0.95 confidence level) as a function of void diameter (fig. 12(b)). Note 
that each point on the plot corresponds to the largest void diameter contained 
In the Interval over which probability of detection was calculated. Hence, 
although the size of the statistical sample was Increased, the calculated prob- 
ability of detection was conservatively Influenced by the Inspection results 
from smaller voids. The optlmlzed-probablllty method demands great computa- 
tional effort, but Its use Is necessary to overcome the problem of Insufficient 
sample size. 



APPENDIX B 


FORTRAN COMPUTER PROGRAM LISTING 

A Fortran computer program that calculates and plots NDE flaw Inspection 
reliability Is listed In this appendix. The program was written for a Digital 
POP 11/45 minicomputer Interfaced with a Grlnnell 274 Image processor. Its 
task file requires 64 kilobytes (125 blocks) of auxiliary memory. Initially, 
flaw Inspection data (the number and size of flaws examined and detected (lee 
table II) Is entered Into the computer and stored In files (or retrieved If 
already stored). This data Is arranged In Intervals by using the options of 
equal flaw size Interval, overlapping flaw size Interval, or optimized prob- 
ability methods (ref. 17). The user selects the Interval sizes and method of 
data grouping. During the grouping, a record of the number and size of flaws 
examined and detected per Interval Is stored. The probability of detection 
pH Is then calculated (ref. 20) over the range of flaw size data using 
any preselected confidence level. Two types of plots displaying quantitative 
flaw detectability results are generated on a video monitor using the Image 
processor software routines. One plot shows the number of flaws examined and 
detected as a function of flaw size or flaw size/part thickness (see figs. 13 
to 20, appendix C). The second plot shows the probability of detection p 4 
(at the selected confidence level) as a function of flaw size or flaw 
size/part thickness (see figs. 8 to 10). 



ORIGINAL PAGa. 

OF. POOR QUALITY 


CALCULATION AND PLOTTING OF NOE FLAW 

INSPECTION RELIABILITY 

Don J. Roth 
August, 1904 


C 

c 


c 


c 

c 

c 


INTEGER Zt YfXf JR (200) fLXfLR(IOO) r LEM 350 > f KB <350) fXYDATA<500> 

1 fXYDAT(SOO) fLLBN< 100) fLLLBN< 100) fLLJZ3<100) 

BYTE ICHAR<2> f IFLPAT<20> 

DIMENSION 1 10 < 300 > » KKQ< 300 ) f AR1 1 < 300 ) t LLR ( 300) f 1X0(300 ) fKXQ<300> 
CHARACTERS J1 f J2*25f J3*2f J4*4f J5*26f J11*5f J22*13f J 33*10 
1 r J44*20 f J6*20f J 50*l f J51*2f J52*2f J53*2f JS4*2f J55*2» JS6#2 
J, f.J57*2f J58#2f ,J59*2f J60#3f J61*3f J62*3f J63*3f J64#3f J65*3 
t 1 J66*3f J67*3» J68*3f J69#3f J70*3f J1 11#17f J1 12*21 f J1 13*14 
If J1 14*20 fJ115*21fJ7*28f J 71*l fJ72*2fJ73*2fJ74*3 
1 f J75*3f J76*3f J77*3f J78*3f J79*3f JB0*3f J81*3f JB2*3f J83*3 
1 9 J84*3 r J85*3 f Je6*3 f J87*3 f J88*3 f J09*3 f J90*3 f J91 *3 f J500* 1 1 
1 F J501*l F J502*2 f J503*2f J900*16f J504*5f J505*5 


Jit**' DATE l ' 

J22=' SAMPLES I.D.J' 

J33~ ' # SCANNED)' 

J44='AV£« THICKNESS (MM) l' 

J5 = ' X NOITCETED FO YTILIBADORP ' 

J6 = ' FLAW SIZE (MICRONS)' 

J7= ' THICKNESS SENSITIVITY (*100)' 

J1 1 l»' CONFIDENCE LEVELS' 

Jl 12« 'DATA GROUPING METHOD I ' 

1113= ' EQUAL INTERVAL' 

JU4 = ' OVERLAPPING INTERVAL' 

Jl .1.5= 'OPTIMIZED PROBABILITY' 

JSOOn'SWALF FO * ' 

J504= ' LATOT ' 

J505= ' DNUOF ' 

J9C0 J ' FL • FREQ. DISTR. ' 

DATA J50/'0'/J51/'10'/J52/'20'/J53/'30'/J54/'40'/J55/'50'/ 
1J5A/'60'/J57/'70'/J58/'80'/J59/'90'/J60/'100VJ&1/'110'/ 

1 J62/ '120' / JA3/ '130' / J64/ ' 140 ' / J65/ ' 150 ' /J66/ ' 160 ' / J67/ '170'/ 
1J68/'180'/ J69/ ' 190 ' /J70/ '200' /J71/ ' 0' / J72/ '25'/J73/ '50 ' / J74/' 75 ' / 
lJ75/'100'/J76/'125'/J77/'150'/J78/'175'/J79/'200'/J80/'225'/ 

1 JQ1/ ' 250 ' / J82/ ' 275 ' /J83/ ' 300 ' / J84/ ' 325 ' / JB5/ ' 350 ' /J86/ ' 375 ' / 

1 J87/'400'/J88/ , 425'/J89/ , 450'/J90/'475 , /J91/'500'/J501/'5'/ 

1 J502/' 15'/J503/'2S'/ 


DATA LMX/0/LLhX/0/MMM/1000/N9 9/1000/ 

DATA IFLDAT/'D'f 'M'f'O'f'J'f'C'f'I'f'f'f'S'fM'f'I'f 
1' 'f' 'f' 'f' 'f' 'f' 'f'.'f'D'f'A'f'T'/ 

INITIALIZE GRINNELL HARDWARE & SOFTWARE 

CALL GRINIT 
CALL GRSINI 
CALL GRSRST 
CALL. GRSBFI (0) 

CALL GRFER <4095f4095f0> ! ERASE EVERYTHING ON VIDEO 

CALL GRNIN (IfOfOfO) (CHANNEL 0 TO ALL LOOKUP TABLES 
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no n r.n n n n r; o no 


CALL GRNBY <1,1, 1,1) ! BYPASS ALL LOOKUP TABLES 

TYPE *, ' 

TYPE # , ' ' 

TYPE ENTER TODAYS DATE (I.E. 12/11/54) »' 

READ (5,59Q) J1 
590 FORMAT (AO) 

Kl«LEN(Jl> ! OBTAIN LENGTH OF DATE 

TYPE *, ' 

TYPE ENTER I.D. OF SAMPLES (EX, AF-BIC-4M' 

READ ( 5, 601 ) J2 
<601 FORMAT (A25) 

K2»LEN<J2> (OBTAIN LENGTH OF 1,0. 

TYPE * , ' ' 

TYPE *,' ENTER * OF SAMPLES SCANNED i ' 

READ (5,602) J3 
602 FORMAT <A2) 


<603 


<611 


604 


K3“LEN ( J3 ) (OBTAIN LENGTH OF * * 

TYPE * » ' ‘ 

TYPE *, ' ENTER AVERAGE THICKNESS OF SAMPLES IN MMt ' 
READ (5,603) J4 
FORMAT <A4) 

K4> t LEN ( J4 ) (OBTAIN LENGTH Of THICKNESS * 

TYPE *, ' 


TYPE *, ' 

TYPE *,' ENTER DATA GROUPING METHOD l ' 

TYPE *,' 0 - EQUAL INTERVAL (OX OVERLAP) POD' 

TYPE 1 - OVERLAPPING INTERVAL POD' 

TYPE *,' 2 - OPTIMIZED PROBABILITY POD' 

TYPE #»' 3 - FLAW FREQUENCY DISTRIBUTION ONLY' 

READ (5,604) L99 
FORMAT (ID 

TYPE *, ' ' 

IF (L99.GT.3) GOTO 611 

IF (L99.EQ.3) GOTO 2153 I* OF FLAWS PLOT 

IF (MMM.NE.1000) GOTO 2178 I PROCEED W/ANOTHER GROUPING METHOD 


629 TYPE *»' ENTER CONFIDENCE LEVEL (X) DESIRED FOR P.O.D.' 
TYPE *,' CALCULATION (EX. 95)' 

READ (5,39) IAR3 
39 FORMAT (12) 

TYPE *, ' 


THE FOLLOWING IS A ROUTINE FOR CHANGING INTEGER 
DATA (CONFIDENCE LEVEL, IAR3) TO CHARACTER DATA 
SO IT CAN BE PRINTED OUT ON VIDEO VIA GRINNELL. 

FOR DETAILED EXPLANATION OF HOW ROUTINE ACCOMPLISHES 
THIS, SEE RBOX PROGRAM, 

ITEMP=I AR3 
RIAR3=FL0ATI ( IAR3) 

IE«INT(AL0G10(RIAR3>> 

rCNT=l 

DO 2150 1*F=0,IE 
(G-"IE-IF 

I TEM=I TEMP/10*# IG 
ICHAR< ICNT ) m ITEM+48 
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ITfc'1P=ITEHP-ITEM*10#*IG 
ICNWCNTfl 
21 SC CONTINUE 

ICNT*ICNT-1 . 

IF (MMM.NE. 1000) GOTO 2178 


QtNAL PALic * 
° F POOR QUALITY 


! PROCEED W/NEW CONFIDENCE LEVEL 


L 

C 

2153 TYPE ENTER I' 

TYPE 0 - TO ANALYZE RADIOGRAPHY DATA' 

TYPE ANY OTHER * - TO ANALYZE SLAM DATA' 

READ (5 . 21 57) L32 
2157 FORMAT (ID 

TYPE *r ' 

C 


C 

c 

c 

C 2178 


2179 


2985 


2987 


2180 

2181 


C 

2281 


686 


687 


786 

C 

688 


C 


THE FOLLOWING ARE GRINNELL MEMORY FILL ROUTINES 
FOR THE P.O.D. VS. FLAW SIZE PLOT TO BE PRINTED 
OUT ON VIDEO 

CALL GRFER (4095r40?5>0) 

CALL GRNIV <1» 1,1,1) ! INVERT GRAPH COLORS 

CALL GRFCD ( 1 r 255 1 0 . 0 . J22 1 10 r 470* 6 * Or 13 , OrO» 0 ) (SAMPLE I.D. 
CALL GRFCDS < 1 t 255 r J2 . 85 , 470 . 6 » 0 r K2 ) I SAMPLE I.D. 

CALL GRFCDS ( 1 , 255 f J33 » 10 » 460 . 6 , 0 r 10 ) I* OF SAMPLES 

CALL GRFCDS ( 1 » 255 . J3 » 80 . 460 7 6 r 0 f K3 ) !* OF SAMPLES 

CALL GRFCDS ( 1.235, J44. 45+K2*6.470. 6,0.20) (THICKNESS 

CALL GRFCDS ( 1 . 255 . J4 » 170+K2#6 . 470 . 6 . 0 .K4 > (THICKNESS 

IF (N99.NE.1000) GOTO 352 (RETURN TO DATA RETRIEVAL ROUTINE 

CALL GRFCDS < 1 . 255 . J1 1 . 202+K2*6+K4*6 . 470 . 6 . 0. 5 ) I DATE 

CALL GRFCDS ( 1 , 255 • JD 242+K2#6+K4*6 . 470 . 6 . 0 t K1 ) (DATE 

IF (MMM.EQ.3.0R.L99.EQ.3) GOTO 2985 It OF FLAWS ON Y-AXIS 

CALL GRFCDS ( 1 , 255 . J5 . 10 . 100 » 0 » 9 t 26 ) (P.O.D 

GOTO 2987 

CALL GRFCDS ( 1. 255 1 J500 . 10 .1 50 . 0 . 9, 1 1 ) !♦ OF FLAWS 

CALL GRFCDS ( 1 . 100. J504.5r2E0.0f 9»5) I TOTAL 

CALL GRFARS ( 1 . 100 t 6 t 310 t 3 . 20 ) (TOTAL-DISPLAY COLOR 

CALL GRFCDS ( 1 » 255 1 J505 . 15 » 250 . 0 .9 . 5 > I FOUND 

CALL GRFARS ( 1 *255. 16 >310 .3.20) (FOUND-DISPLAY COLOR 

IF (L32.EG.0) GOTO 2180 I RADIOGRAPHY 

CALL GRFCDS < 1 . 255 . J6 . 200 . 10 . 6 » 0 . 20) (FLAW SIZE 

GOTO 2181 

CALL GRFCDS < 1 » 255 » J7» 200 t 10 . 6 . 0 » 28 > (THICKNESS SENSITIVITY 
CALL GRFCDS ( 1. 255 . J1 1 1 . 100. 460 t 6, 0 . 17 > (CONFIDENCE LEVEL 

IF (L99.ER.3) GOTO 2281 

CALL GRFCDS ( 1 .255. ICHAR.206 , 460.6 » 0/ ICNT> (CONFIDENCE LEVEL 

CALL GRFCDS ( 1 .255. J1 12 r 240» 460 » 6.0.21 > (DATA GR. METHOD 
IF (L99.EQ.0) GOTO 686 (EQ. INT. 

IF (L99.EG.1) GOTO 687 ! OVERL . INT. 

IF (L99.E0.3) GOTO 786 I# OF FLAWS PLOT 

CALL GRFCDS ( 1 » 255 ► J1 15 . 380. 460 . 6 . 0 . 21 ) (OPT. PROB. 

GOTO 688 

CALL GRFCDS ( 1 . 255 * J1 13 . 380 f 460 * 6 . Or 14 ) ! EQ . INT. 

GOTO 688 

CALL GRFCDS ( 1 . 255 r J1 1 4 . 380 . 460. 6 .0.20 ) (OVERL. INT. 

GOTO 688 

CALL GRFCDS ( 1 . 255 . J?00 . 380 » 460. 6 . 0 . 16 ) (# OF FLAWS PL. 

CALL GRFVC ( 1 .255.0.0.50.50.450,50) (X-AXIS 
CALL GRFVCS ( 1 , 255 , 50, 50 , 50 , 450) (Y-AXIS 
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i* OF FLAWS PLOT 


IF ( MMM . EQ . 3 . OR * L99 • EQ . 3 > GOTO 68? 

GOTO 6161 I PROBABILITY PLOT 
689 HO 6121 119=50, 150,16 I* OF FLAWS PLOT Y HASH MARKS 

CALL GRFVCS (1 255,45 ,119,50,119) 

6121 CONTINUE 

IF (L32.EQ.0) GOTO 6026 I RADIO. F\, FREQ. X-AXIX HASH 
GOTO 6311 I X-AXIS HASH MARKS FLR SLAM FL . FR. PL. 

C 

6161 DO 625 J 19=50 , 450 , 8 

CALL GRFVCS <1.255.45.119.50.119) (Y-AXIS (SMALL) HASH MARKS 

625 CONTINUE 

DO 6025 122=50,450,40 

CALL GRFVCS (1,255,40,122,50,122) IY-AXIS (LARGE) HASH MARKS 

6025 CONTINUE 

IF (L32.EQ.0) GOTO 6026 IRADIUG. 

6311 INCR1=10 I FLAW SIZE HASH MARKS 

1NGR2=20 (DITTO 

GOTO 6131 

6026 INCR1=9 (THICKNESS SENSITIVITY HASH MARKS/RADIO. 

INCR2=40 I DITTO 

6131 DO 626 120=50 , 450, INCR1 

CALL GRFVCS (1,255,120,45,120,50) (X-AXIS (SMALL) HASH MARKS 

626 CONTINUE 

DO 627 121=50,450, INCR2 

CALL GRFVCS (1,255,121,40,121,50) (X-AXIS (LARGE) HASH MARKS 

627 CONTINUE 

IF ( MMM . EQ . 3 . OR , L99 , EQ , 3 ) GOTO 3864 !♦ OF FLAWS PLOT 

C 

C 

CALL GRFCDS ( 1 , 255 , J50 , 35 , 50 , 6 , 0 » 1 ) (Y-AXIS *'S 

CALL GRFCDS ( 1 , 255 , J51 , 30 , 90 , 6 , 0 , 2 > 

CALL GRFCDS ( 1 , 255 , J52 , 30 , 1 30 , 6 , 0 , 2 ) 

CALL GRFCDS ( 1 , 255 , JS3 , 30 , 170 , 6 , 0 , 2 ) 

CALL GRFCDS ( 1 ,255, J54 , 30, 210, 6 ,0, 2 ) 

CALL GRFCDS ( 1 , 255 , J55 , 30 , 250 , 6 , 0. 2 ) 

CALL GRFCDS ( 1 , 255 , J56 , 30 , 290 , 6 , 0» 2 ) 

CALL GRFCDS ( I , 255 , 057 , 30 , 330 , 6 , 0, 2 ) 

CALL GRFCDS ( 1 , 255 , J58 , 30 , 370 , 6 ,0 ,2) 

CALL GRFCDS ( 1 , 255 , J59 , 30 , 410 , 6,0 . 2 ) 

CALL GRFCDS ( 1 , 255 , J60, 25 , 450 , 6,0 , 3) 

CALL (3RSDFD 
GOTO 3865 

3864 CALL GRFCDS . i , 255 , J50 , 35 , 50 , 6 , 0 , 1 ) !* OF FLAWS Y-AXIS #'S 

CALL GRFCDS ( 1 , 255 , J501 , 35 , 1 30 , 6 , 0 , 1 ) 

CALL GRFCDS ( 1 , 255 , J51 , 30, 210 , 6,0 , 2) 

CALL GRFCDS ( 1 , 255, JS02, 30, 290,6 ,0, 2) 

CALL GRFCDS <1 ,255, J52, 30, 370, 6, 0,2) 

CALL GRFCDS ( 1 ,255, J503, 30 ,450,6 ,0,2 ) 

CALL GRSBFD 

3865 IF (L32.EQ.0) GOTO 628 (RADIOGRAPHY 
C 

CALL GRFCDS (1,255,050,50,30,6,0,1) (X-AXIS *'S 
CALL GRFCDS (1,255,051,65,30,6,0,2) 

CALL GRFCDS (1,255,052,85,30,6,0,2) 

CALL GRFCDS ( 1 , 255 ,• 053 , 105 , 30 , 6 , 0 , 2 ) 

CALL GRFCDS ( 1 , 255 , 054 i 1 25 , 30, 6 , 0 , 2 ) 

CALL GRFCDS (1,255,055,145,30,6,0,2) 

CALL GRFCDS ( 1 ,255, 056 , 165 ,30,6 , 0,2) 

CALL GRFCDS ( 1 , 255 , 057 , 185, 30, 6 , 0 , 2 ) 

CALL GRFCDS ( 1 , 255 , 058 , 205 , 30 , 6 , 0 , 2 ) 

CALL GRFCDS (1,255,059,225,30,6,0,2) 
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CALL GRFCDS < 1 p25S p J60 p 242, 30 r At Or 3 ) 

CALL GRFCDS ( 1 p255f J61 f262p30p6f0p3) 

CALL GRFCDS < 1 1 255r J62? 2B2 p 30 p 6 p Or 3 ) 

CALL GRFCDS < 1 r 255 r J63 t 302 f 30 r 6»0r 3 > 

CALL GRFCDS ( 1 f 255 r J64 r 322 r 30 p 6 tO t 3 > 

CALL GRFCDS ( 1 r25Sr J65r342r30r6r0f3> 

CALL GRFCDS ( 1 r 255 r J66 r 342 r 30 f 6 ? Or 3) 

CALL GRFCDS ( 1 r 255 r J67 r 382 r 30 . 6 r Of 3 > 

CALL GRFCDS ( 1 r 255 r J6ep402r30r6«0r3> 

CALL GRFCDS ( 1 r 255 » J69 » 422 p 30 r 6r0r 3 ) 

CALL GRFCDS < 1 r 255 r J70 r 442 r 30 p 6 r 0 p 3 > 

CALL GR8BFD 
GOTO 59 
C 

628 CALL GRFCDS < 1 r 255 t J71 t 50 » 30 r 6 p 0 p 1 ) (RADIOGRAPHY X-AXIS *'S 

CALL GRFCDS < 1 r255r< J73r85r30r 6 f0p2) 

CALL GRFCDS ( 1 r 255 r J75 r 122 r 30 r 6 rO r 3 > 

CALL GRFCDS < 1 r255, J77 p 162p30p6p0p3) 

CALL GRFCDS < 1 p255p J79p202»30p6pOp3) 

CALL GRFCDS <1 p255p J 81 p 242 p 30 p 6 p 0 p 3 ) 

CALL GRFCDS (1 . 255 p J83 p 282 p 30 p 6 r 0 p 3 > 

CALL GRFCDS < 1 p 255 t J85 p 322 p 30 p 6 t 0 p 3 ) 

CALL GRFCDS < 1 p 255 p J87 » 362 t 30 p 6 f 0 r 3 ) 

CALL GRFCDS ( 1 p 255s J89 p 402 p 30 p 6p0 p 3 > 

CALL GRFCDS < 1 p 255 p J 91 p 442 p 30 p 6 p 0 p 3) 

CALL GRSBFD 

t 

59 IF ( MMM . EQ . 0 . AND . L99 • EG > 0 > GOTO 219 ! EQ . INT. PATH 

IF (MMM.EQ. 0.AND.L99.EQ. 1 ) GOTO 219 (QVERL. INT. PATH 
IF (MMM.EQ. 1 . AND.L99.EG.2) GOTO 6001 IOPT. PROS. PATH 

IF (MMM.EQ.2. AND.L99.EQ.2) GOTO 2005 (OPT. PROD. PATH 

IF (MMM . EG. 0 . AND . L99. EG. 2 ) GOTO 2005 (OPT. PROD. PATH 

IF (MMM.EG.1.AND.L99.EG.0) GOTO 500 (EG. INT. PATH 

IF (MMM.EG.l .AND.L99.EG, 1 > GOTO 500 ! OVERL . INT. PATH 

IF < MMM . EG . 2 . AND . L99 . EQ . 0 ) GOTO 219 (EG. INT. PATH 
IF ( MMM . EG . 2 . AND . L99 . EQ . 1 ) GOTO 219 (OVERL. INT. PATH 
IF (MMM. EG. 3) GOTO 6028 I# OF 1 LAWS PLOT 

INITIALIZE FLAN SIZE DATA ARRAYS 

DO 61 MB“1p600 
XYDATA(MEO-0 
XYDAT (MB)=0 
61 CONTINUE 

DO 8 MD a l p 300 
LD(MB)-0 
KB ' MB ) “0 
8 CONTINUE 


TYPE *>' 

TYPE *r' ENTERS ' 

TYPE *r' 0 - TO ENTER NEW FLAW SIZE DATA' 

TYPE 1 - TO RETRIEVE STORED FLAW SIZE DATA' 

READ ( 5 p 343 ) N99 
343 FORMAT (ID 

TYPE *» ' 

IF (N99.EG.0) GOTO 378 (ENTER NEW FLAW SIZE DATA 


THE FOLLOWING ROUTINE RETRIEVES FLAW SIZE DATA FROM DISK i 
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c (FOR DETAILED DOCUMENTATION OF STORING AND RETRIEVING DATA * 

C SEE RBOX PROGRAM) 

C 

TYPE 

TYPE TYPE IN THE 4-CHARACTER NAME OF DATA FILE YOU WISH" 

TYPE TO RETRIEVE f (IF YOU INPUT 0CZER03r RETRIEVAL WILL BE' 
TYPE ABORTED AND PROGRAM WILL PROCEED W/NEW DATA ENTRY' 
ACCEPT 350 » IFLDAT (11) , IFLDAT ( 12) , IFLDAT ( 13 ) » IFLDAT ( 14 ) » 

1 IFLDAT ( 15 ) f IFLDAT (14) 

350 FORMAT (4A1) 

TYPE *,' ' 

IF (IFLDAT (11) .EQ, '0') GOTO 378 

OPEN ( UNI T«4 , NAME® I FLOAT , STATUS® ' OLD ' , FORM® ' UNFORMATTED ' , 

1 BLOCKSIZE<=1024) 

r, 

C RETRIEVE FLAW DATA FORM DISK 

C 

READ ( UNI T=4 ) J4 , I B , LB , KB , NAN , LLBN » LL JZ3 , LLLBN 
CLOSE (4) 

GOTO 2179 I TO PRINT THICKNESS ON VIDEO 

352 TYPE ' 

TYPE ENTERS' 

TYPE *r' 0 - WRITE 0U T FLAW SIZE FILE DATA FILE TO TERMINAL' 

TYPE *, ' ANY OTHER * - DO NOT WRITE FILE TO TERMINAL' 

READ (5,3052) IJIJ 
3052 FORMAT (II) 

IP (lJIJ.NE.O) GOTO 444 
WRITE (5,355) IFLDAT, IB 
355 FORMAT (l.OX , 20A1 , 5X, 13 , IX, 'ENTRIES ' ) 

TYPE *, ' ' 

DO 353 N84«i»IH 

WRITE (5,354) N84,LB(N84) ,KD(N86) 

354 FORMAT ( 10X , 'ENTRY ( ' , 13 , ' ) ' ,2X, 13 , IX, I 1 ) 

353 CONTINUE 


C 

3332 


3333 


C 

c 


3335 

3334 


TYPE 

TYPE 

TYPE 

TYPE 

TYPE 

TYPE 

TYPE 

TYPE 

READ 


/ 


to 

*, 'DO 

*, ' 


*, ' 

*, ' 

*, ' 

*, ' 

#, ' 

(5,3333) 
FORMAT (ID 
IF (NONO.GT.l) 
IF (NONO.EG.l) 


YOU WISH TO CHANGE FLAW SIZE DATA FILE?' 
ENTERS ' 

CHANGE AN ENTRY' 

ADD ENTRIES TO FILE (FIRST CHANGE ' 
LAST ENTRY C00003 OF ORIGINAL FILE' 
TO FIRST NEW ENTRY) ' 

OTHER * - NO CHANGE' 

NOND 


0 - 
1 - 


ANY 


GOTO 444 ! NO CHANGE 

GOTO 3884 I ADD ENTRIES 


TYPE 

TYPE ENTER * OF ENTRY YOU WISH TO CHANGE (13 FORMAT)' 

READ (5,3335) JIB 
FORMAT (13) 

WRITE (5,3334) JIB,LB( JIB) ,KB( JIB) 

FORMAT (IX, 'FLAW SIZE IX, 13 IX, 13, IX, II ) 

TYPE *, ' 

Type tr ' ENTERS ' 

TYPE *,' 0000 - TO DELETE FLAW SIZE ENTRY' 

TYPE *, ' 13,11 # - REPLACEMENT FLAW SIZE ENTRY' 
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nn 


READ (5,3337) LB< JIB) »KB< JIB) 
3337 FORMAT <I3,IL) 

GOTO 3993 


3884 TYPE ENTER # OF NEW ENTRIES YOU WISH TO ADD TO FILE' 

TYPE (NOTE! OOOO MUST BE INCLUDED AS FINAL ENTRY)' 

READ (5,3885) IBP 

3085 FORMAT (13) 

DO 3886 NET*1 , IBP 
WRITE (5,3887) NET 

3887 FORMAT (SX, 'ENTRY • ' , IX, 13 , IX , ' * ( 13, I 1 FORMAT)') 

READ (5,3880) LB ( IB+NET ) ,KB ( IB+NET ) 

3888 FORMAT (13,11) 

3086 CONTINUE 

I B ■»IB + IBP 
TYPE ' 

2993 OPEN ( UNI T*4 , NAME* IFLDAT , STATUS* ' UNKNOWN ' , FORM* ' UNFORMATTED ' , 

1 BL0CKSIZE=1024> ! PUT CHANGE IN FILE 

WRITE (UNIT=4> J4, IB,LB,KB,NAN,LLBN,LLJZ3»LLLBN ! WRITE CHANGE 

CLOSE (4) 

C 

OPEN ( UN I T*4 , NAME* IFLDAT , STATUS* ' UNKNOWN ' , FORM* ' UNFORMATTED ' , 

1 BLOCKS I ZE = 1024 > (TAKE CHANGED PILE AND RECALCULATE 

C 

C 

(, PERFORM OPERATIONS ON FILE DATA AS WAS DONE TO 

C INITIALLLY-ENTERED DATA >ivT nov*i 


READ <UNIT*4) J4 , IB , LB , KB , NAN , LLBN , LL JZ3, LLLBN 

CLOSF (4) 

LZZ a G 

NAN~0 

JZ3“0 

LB ( 0 ) »LB ( 1 ) 

DO 8490 JZF=1 , ID 

IF (LB< JZF) .EQ.LB( JZF-1) ) GOTO 8492 

NAN*NAN+1 

LLBN ( NAN >»LZZ 

LLJZ3 ( NAN ) JZ3 

LLLBN ( NAN )*LB( JZF-1) 

IF (JZF, ED . IB) GOTO 8490 (END RECALCULATION OF FILE DATA 

LZZ = 1 

IF (KB(JZF).NE.l) GOTO 8495 
JZ3-1 
GOTO 8490 
8495 JZ3«0 

GOTO 8490 

8492 IF (KB (JZF) ,NE.l) GOTO 8493 
JZ3-JZ3+1 
0493 LZZ-LZZil 
8490 CONTINUE 
C 

C PUT RECALCULATED LLBN, LLJZ3, LLLBN ARRAYS IN FILE 

C 

, OPEN ( UNIT“4 , NAME=IFLDAT , STATUS* ' UNKNOWN' , FORM* 'UNFORMATTED ' , 
1 BLOCKS I ZE* 1024 ) 

WRITE ( UNIT=4 ) J4 , IB, LB , KB , NAN, LLBN,LLJZ3,LLLBN 
CLOSE (4) 

OPEN <UNIT=4,NAME*IFLDAT * STATUS* ' UNKNOWN ' , FORM* ' UNFORMATTED ' » 
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3334 


C 

C 

0 

c 

c 


378 


379 


400 

401 


402 


403 

403 


404 

407 


C 

r; 

454 


458 


C 

c 

c 


1 BLOCKS I ZE» 1024) I READ OUT NEW FILE 

READ <UNIT«4) J4 * IB * LB *KB * NAN* LLBN* LLJZ3 * I.LLBN 
CLOSE (4) 

IF ( NONQ • EO . 0 . OR . NONO * EQ . 1 ) GOTO 3332 

TYPE ** ' 

GOTO 464 I SKIP NEW FLAW BATA ENTRY 


THE FOLLOWING ROUTINE ALLOWS THE USER TO ENTER 
FLAW SIZE BATA OBTAINED ON METALLOGRAPH 

IB=0 ! INITIALIZE TOTAL FLAW ♦ COUNTER 
NAN“0 (INITIALIZE * OF SPECIFIC FLAW SIZES COUNTER 
LZZ“0 (INITIALIZE TOTAL FL.AWS/SP, SIZE COUNTER 
JZ3«0 (INITIALIZE FLAWS FOUND/SP. SIZE FLAW COUNTER 
TYPE #* ' 

TYPE #*' ENTER FLAW BATA FROM METALLOGRAPH' 

TYPE **' AND SLAM/RADIOGRAPHY <13*11 FORMAT)' 

TYPE *r ' ENTER * ' 

TYPE **' 1) FLAW SIZE FROM METALLOGRAPH IN ASC. OR DESC. 

TYPE #* ' 2) 1 - IF FLAW FOUND ON SLAM' 

TYPE **' 2 - IF FLAW NOT FOUND ON SLAM' 

TYPE * * ' ( ENTER 0<ZER0) AFTER LAST FLAW IS ENTERED)' 

TYPE ** ' ' 

READ <5*401) JB*JU 
FORMAT (13*11) 

LB(0)“LB(1> 

IB-IB+1 (IBM OF FLAWS ENTERED 

LB( IB)=JB (STORE FLAW SIZES IN ARRAY 
«B(IB)=JU (STORE SLAM RESULT IN ARRAY 
IF <LB(IEO .EG.LB(IB-l) ) GOTO 402 
GOTO 403 

IF (KB (IB) .NE.l) GOTO 405 

JZ3= JZ3+1 (KEEP TRACK OF # OF FLAWS FOUND/FL. SIZE 
LZZ l -LZZ+l I KEEP TRACK OF # OF SAME FLAWS 
GOTO 400 
NAN M NAN+ 1 

LLBN ( NAN ) =LZZ t* OF FLAWS OF PARTICULAR SIZE 
I,.L JZ3 < NAN ) = JZ3 I# OF FLAWS FOUND OF PARTIC. SIZE 
LLLDN( NAN) “LB < IB-1 ) (PARTICULAR SIZE FLAW/TH. SENS. 

LZZ=*1 

IF ( KB ( IB > . NE • 1 ) GOTO 406 

JZ3=1 

GOTO 407 

JZ3=0 

IF (JB.EO.O) GOTO 454 
GOTO 400 


TYPE ** ' 

TYPE #* ' ' 

TYPE **' ENTER*. ' 

TYPE **' 0 - TO FIRST STORE DATA IN FILE AND THEN GROUP' 
TYPE **' ANY OTHER * - PROCEED DIRECTLY TO DATA GROUPING' 
READ < 5*450) LM5 
FORMAT <I1) 

TYPE ** ' 

IF (LM5.NE.0) GOTO 464 IDO NOT STORE DATA 

THE FOLLOWING ROUTINE STORES FLAW SIZE DATA ON DISK 
(FOR DETAILED DOCUMENTATION OF STORING AND RETRIEVING 


ORDER' 
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c 

c 


DATA, SEE RBOX PROGRAM) 


TYPE *f ' 

TYPE TYPE IN A ^-CHARACTER NAME OF DATA FILE YOU WISH' 

TYPE TO STORE f (IF YOU INPUT 0CZERQ3 1 STORAGE WILL BE' 

TYPE *>' ABORTED AND PROGRAM WILL PROCEED W/ DATA GR,)' 

TYPE * r ' 

ACCEPT 459 f IFLDAT (11) , IFLDAT ( 12) » IFLDAT ( 13) f IFLDAT( 14) t 
1 IFLDAT < 15 ) » IFLDAT ( 16 ) 

459 FORMAT <6A1> 

TYPE %, ' ' 

IF ( IFLDAT ( 1 1 ) .EQ. ' 0 ' > GOTO 464 

OPEN < UN I T“4 » NAME* IFLDAT t STATUS" ' UNKNOWN' rFORM-'UNFORMATTED' r 
1 BL0CKSIZE=1024> 

C 

C STORE FLAW SIZE DATA ON DISK 

C 

WRITE (UNIT=4) J4 * IB f LB r KB » NAN rLLBN f LLJZ3 oLLLBN 

CLOSE (4) 

TYPE *t ' ' 

WRITE (5» 437) IFLDAT 

437 FORMAT < 5X , 20A1 > 1 X » ' HAS BEEN STORED') 

TYPE *r ' ' 

C 

C 

464 IF (MMM.EQ.3.0R.L99.EQ.3) GOTO <4020 !♦ OF FLAWS PLOT 

IF (L99.EQ.0 .OR. L99.EQ.1) GOTO 219 ! OVERLAPPING/EQUAL INT. 

GOTO 2005 I OPT. PROB. ROUTINE 

6028 IF (L32.EQ.0) GOTO 6040 (RADIOGRAPHY 
C 

C 

C 

C FLAW FREQUENCY ROUTINE 

C 

C 

DO 6030 1.57=1 f NAN I* OF FLAWS PLOT/SLAM DATA 

IF (LLBN(L57) .LE.27) GOTO 6029 

LLBN(L57>=27 

IF (LLJZ3(L57> .LE.27) GOTO 6029 
LLJ"3(L57)=27 

6029 CALL GRFUCS ( 1 r 100r50+LLLBN(L57)#2>50+LLBN(L57)*16* 

150+.-LLBN (L57)#2»50) (TOTAL FLAWS OF PARTIC. SZ. 

CALL GRFVCS ( 1 r255»50+LLLBN(L57)#2fS0+LLJZ3(L57>*16» 
1504-LLLBN(LS7)#2»50> (FLAWS FOUND 

6030 CONTINUE 

GOTO 8020 (END PROGRAM 

6048 DO 6050 L57=l fNAN !* OF FLAWS PLOT/RADIOGRAPHY DATA 
IF (LLBN(L57) , LE.27) GOTO 6049 

LLBN ( L57 ) =27 

IF (LLJZ3(L57) .LE.27) GOTO 6049 
LLJZ3(L57>=27 

6049 CALL CRFUCs" ( 1 , 100» 50+NINT(LLLBN(L57)» .0) »50+LLBN(L57>*16» 

1S0+NINT(LLLBN(L57>*,8) ,50) !* OF FLAWS OF PARTIC. SZ, 

CALL GRFVCS ( 1 ,255r50+NINT(LLLBN(L57>*.8> r50+LLJZ3(L57>*16» 
150+NINT(LLLBN(L57>*.8) »50) !# OF FLAWS FOUND 

60.50 CONTINUE 

CALL GRSBFD 

GOTO 8020 (END PROGRAM 
C 
C 


21 



n o > - o on 


OPTIMIZED PROBABILITY METHOD FOR GROUPING DATA 

2003 TYPE *, ' ' 

DO 71 M0P«1 , 300 UNITIALIZE ALL ARRAYS 
LLR ( MOP ) ®0 

iiq(Mop>»o 

IXQ(M0P><=0 
KKQ<M0P>«0 
KXQ < MOP ) °0 
ARU < MOP ) “0 
71 CONTINUE 

TYPE ENTER t ' 

TYPE 0 - TO WRITE ALL CALC. PROB. TO TERMINAL' 

TYPE ANY OTHER # - TO OMIT WRITING CALC. PROB,' 

READ ( 5 , 2007 ) MOM 
2007 FORMAT (ID 

TYPE #, ' ' 

TYPE *,' ENTER i' 

TYPE 0 - PLOT FLAW SIZE RANGE BARS' 

TYPE ANY OTHER * - DO NOT PLOT FLAW SIZE RANGE BARS' 

READ (5,907) MAZE 
987 FORMAT (ID 
TYPE ' 

TYPE *, ' ' 

TYPE *, ' ' 

TYPE ENTER THE # OF BOUNDARY FLAW SIZES TO BE USED!' 

READ ( 5,2010 ) IIX 
2010 FORMAT (12) 

TYPE *, ' ' 

TYPE IMPORTANT)!' 

TYPE #»' NOTE! PLEASE TYPE BOUNDARY FLAW SIZES IN' 

TYPE *,' ASCENDING ORDER' 

TYPE ' ' 

DO 3050 JJX=1,IIX 
WRITE (5,4000) JJX 
4000 FORMAT (' BOUNDARY FLAW SIZE t ' , 1 3 , ' = ( ENTER SIZE IN MICRONS) ') 
READ (5,5555) LLX 
5555 FORMAT (13) 

LLR ( JJX ) =LLX (STORE BOUNDARY FLAW SIZES IN ARRAY 

3050 CONTINUE 
TYPE 

TYPE ' 

THIS ROUTINE PLACES FLAWS IN PROPER INTERVALS 


IIQ( 1) =0 (INITIALIZE ♦ OF FLAWS/INTERVAL COUNTER 
IXQ ( 1 ) =0 (SAME AS ABOVE FOR IXQ 

KKQ( 1 ) =0 (INITIALIZE « OF FLAWS FOUND/INTERVAL COUNTER 
KXQ( 1 ) =0 (SAME AS ABOVE FOR KXQ 
JJQ=1 (SET INTERVAL COUNTER INITIALLY TO 1 

3057 DO 5900 NNB=1,IB (FOR EACH FLAW 

IF (LB(NNB) .GE. LLR (JJQ) .AND.LB(NNB) »LE.LLR( JJQ+1 ) ) GOTO 5050 
GOTO 5900 

5050 Iia(JJQ)=IIG(JJQ)+l (KEEP TRACK OF * OF FLAWS/INTERVAL 

IXQ ( J JQ ) =1X0 ( JJQ > +1 (SAME AS ABOVE FOR IXQ ARRAY 

IF (KB(NNB) .EG. 2) GOTO 5900 (FLAW NOT SEEN ON SLAM 
KKQ ( J JQ ) =KFQ ( J JQ ) +1 (KEEP TRACK OF * OF FLAWS SEEN/INTERVAL 

KXQ ( JJR ) -”KXQ( JJQ) + 1 (SAME AS ABOVE FOR KXQ ARRAY 

5900 CONTINUE 

WRITE (5,5910) KKQ ( J JQ ) , I IQ ( JJQ ) ,LLR( J JQ ) , LLR ( JJQ+1) 
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ORIGINAL PAGE- is 
of POOR QUALITY 


5910 FORMAT ( 1 X, 13 , IX * 'OUT OF ' » 13 , 1 X , 'FLAWS BETWEEN '.IX, 14, IX 
1 , ' AND' , 14 » 1 X , 'MICRONS WERE FOUND ON SLAM') 

TYPE 

JwlQraJJG+1 (INCREMENT INTERVAL COUNTER 

IF ( JJQ • EO i I IX) GOTO 6001 (PROCEED TO PROB. CALCS, 

GOTO 3057 (TRY ALL FLAWS AGAIN W/NEW INTERVAL 

C 

C THIS ROUTINE CALCULATES PROBABILITIES <OPT. PROP.) 

C 

r 

" 4001 MMQ“I IX-1 I SET MMQ COUNTER INITIALLY TO • OF INTERVALS 

C 

6100 JJO^MMQ (SET JJQ TO MMQ 

DO 6200 LED«1 , IIX-1 (FOR EACH INTERVAL 

I IQ ( LED ) “IXCHLED > 'RESET IIQ VALUE FOR CORRECT CUM, PROP. CALC, 
KKQ < LED ) “=KXQ ( LED ) (RESET KKQ VALUE FOR CORRECT CUM, PROB, CALC, 
6200 CONTINUE 
C 

6500 AR1 k FL0ATI(IIQ(JJQ)) (CHANGE IIQ TO REAL 
AR2>=FL0ATKKKQ< JJQ) > (CHANGE KKQ TO REAL 

AR3«FL0ATI(IAR3>/100, (CHANGE CONF, LEVEL TO REAL/DECIMAL 
CALL DIN< AR1 » AR2 , AR3 , AR10 > (SUBROUTINE FOR PROB. CALC, 

ARll < JJQ)«AR10 
IF (MOM.NE.O) GOTO 6502 
WRITE (5.6501) JJQ » AR1 1 ( JJQ ) 

6501 FORMAT (lOX.'CALC, PROB. < ' . 12 . ' )= 'F5 , 3 ) 

6502 IF ( JJQ , EG , 1 > GOTO 7000 (FINISHED W/THIS INTERVAL 

JJQ 33 JJQ- 1 (DECREMENT J IQ BY 1 

I ICK JJQ ) =1 IQ ( JJQ+1 > +1 IQ ( JJQ) (CUMULATIVE TOTAL FLAWS 

KKQ< JJQ> 33 KKQ( JJQ+1)+KKQ<JJQ) (CUMULATIVE FLAWS FOUND 
GOTO 6500 


' ROUTINE TO GET MAX PROBABILITY FROM ABOVE CALC. PROBABILITIES 

C 

7000 PMAX»0, 

DO 7020 JKQXbI'MMQ I FOR EACH PROBABILITY 
IF (ARll(JKQX) .GT.PMAX) GOTO 7010 I TEST 

GOTO 7020 

7010 PMAX=ARl 1 < JKQX > 

JKQ-JKQX 
7020 CONTINUE 
C 

5920 WRITE (5.5930) IAR3.PMAX 

5930 FORMAT ( IX. 'AT IX . 12. ' X IX. 'CONFIDENCE. THE LOWER-BOUND 
1 PROBABILITY OF DETECTION‘S '» IX. F5. 3) 

WRITE (5,5940) (LLR(MMQ)+(LLR(MMQ+1 >-LLR(MMQ> ) ) 

5940 FORMAT <20X,'AT FLAW SIZE 33 ' » IX. 13) 

LLMX=LLMX+1 (KEEP TRACK OF * OF PTS. TO BE PLOTTED 

TYPE ' 

IF (L32.EQ.0) GOTO 5943 (PLOT FOR RADIOGR. 

XYDATA ( 2DCLLMX-1 >=50+ (LLR(MMQ)+(LLR(MMQ+1 >-LLR<MMQ) > )*2 (FL.SZ, C. 
XYDAT ( 2DCLLMX-1 ) =50+ (LLR(MMQ) + ( LLR(MMQ+1 )-LLR(MMQ) ) )*2 ITH.PL. 

GOTO 5944 (SKIP RAD, POINTS 

5943 XYDATA ( 2#LLMX-1 ) -50+NINT ( (LLR(MMQ) + (LLR(MMQ+1 ) -LLR(MMQ) ) >4 .0 > (RD. 
XYDAT <2D(LLMX-1 )=50+NINT ( <LLR(MMQ) + (LLR(MMQ+1)-LLR(MMQ) > >*.8> (RD. 

5944 XYDATA (2D<LLMX) =504- <NI NT (PMAXDtlOO. )#4. ) (P.O.D. COOR. 

XYDAT (2DcLLMX)~49+(NINT (PMAXDclOO. )Dt4. ) (THICKEN PLOT 

C 

C PLOT POINTS USING ABOVE COORDINATES (VISIBILITY AID) 

C 
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CALL GRFAR (1 , 254 , 0 , 0 , XYDATA< 24LLMX-1 ) -2»XYDATA< 2*LLMX )-2> 4 ,4 > 

C 

IF (MAZE.NE.O) GOTO 5947 IDO NOT PLOT FL. SZ. RANGE BARS 

C 

C PLOT FLAW SIZE RANGE BARS 

C 

IF (L32.E0.0) GOTO 5946 (PLOT FOR RADIOGR . 

IF ( JKO . EO * 1 * AND • MHO . EO . 1 ) GOTO 4695 

MKG=MMO+l 

GOTO 4696 

4695 MKQ«MHQ+1 

4696 CALL GRFVC ( 1 , 150, 0 , 0, XYDATA<2*LLMX-1 ) -(LLR<HKQ >- 

1LLR( JKO) )*2,XYDATA(2*LLMX> , XYDATA< 24LLMX-1 > , XYDATA < 2*LLMX > > 

C 

CALL GRFVC ( 1 r 150,0, 0, XYDATA ( 24LLMX-1 ) -<LLR (MKO >- 
1LLR ( JKO ) ) *2 , XYDATA ( 2*LLMX ) -2 ,X YDATA< 2*LLMX-1 ) - 
.1 < LLR < MKO ) -LLR ( JKQ) ) *2 , XYDATA < 2*LLMX > +2 ) 

GOTO 5947 I SKIP RADIOGRAPHIC POINTS 

C 

5946 IF ( JKO . EO . 1 . AND • MMQ • EQ . 1 > GOTO 4211 
MKQwMMG+1 

GOTO 4213 
4211 MKQaMMO+1 

4213 CALL GRFVC ( 1 , 150,0 , 0»XYDATA ( 2SLLMX-1 >-NINT< < LLR(MKQ) 

1 -LLR (JKQ) )*.S> ,XYDATA(2#LLMX) ,XYDATA(2*LLMX-1 > , 

1XYDATA<2*LLMX ) ) 

C 

CALL GRFVC < 1 , 150,0,0,XYDATA(2*LLHX» 1 )-NINT< <LLR<MKO) 

1 -LLR (JKQ) )#.8> » XYDATA (2*LLMX) -2, 

1 XYDATA (2*LLMX~1 >-NINT< (LLR(MKQ) -LLR( JKQ) >#.8>, 

!XYDATA(2*LLMX)+2) 

C 

C 

c 

5947 WRITE (5,5949) MKQ, JKQ, LLR < JKQ) , LLR (MKO) 

5949 FORMAT < 10X , 13 , SX , 13 , 5X , 13 , 5X , 13 > 

MMQ°MMG~1 

IF (MMQ.EQ.O) GOTO 8010 
GOTO 6100 

8010 CALL GRFVL ( 1 , 254 » 0 , 0, XYDATA , LLMX) I PLOT OF P.O.D. VS. FLAW SIZE 
CALL GRFVL ( 1 , 254 ,0 ,0 ,XYDAT ,LLMX > (THICKEN PLOT LINE 

CALL GRSBFD 

LLMX-0 

N99 c =1000 

Type; £ f * ENTER J ' 

TYPE 0 - TO GROUP DATA USING EQUAL/OVERL . INT. METHOD" 

TYPE 1 - TO USE DIFFERENT CONFIDENCE LEVEL" 

TYPE 2 - TO USE DIFFERENT INTERVALS' 

TYPE 3 - TO PLOT * OF FLAWS VS. FL.SZ./TH. SENS." 

TYPE 4 - TO END PROGRAM' 

READ (5,8011) MMM 

8011 FORMAT (ID 

TYPE #, " " 

IF (MMM.EQ.O) GOTO 611 
IF (MMM.EQ.l) GOTO 629 
IF (MMM.EQ.2) GOTO 2178 
, IF (MMM.EQ.3) GOTO 2178 
GOTO 8020 
C 

c 

• C END OF OPTIMIZED PROBABILITY METHOD 
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START OF EQUAL/OVERLAPPING INTERVAL METHOD FOR DATA 


21? TYPE *r' ' 

TYPE *•' 

DO 70 HB n l r 100 UNITIALIZE ALL ARRAYS 
LR(MB>»0 
70 CONTINUE 

DO 60 MB«1 » 200 
JR< MB >“0 
60 CONTINUE 

TYPE # r ' ' 

TYPE **' ENTER THE « OF BOUNDARY FLAW SIZES TO BE USED t ' 

TYPE *r' (REMEMBER! BOUNDARY FLAW SIZES MUST BE EQUALLY SPACED)' 
TYPE *t ' (EX * 4Of0Of12Of 160 MICRONS)' 

READ ( 5 f 2 1 0 ) IX 
210 FORMAT (12) 

TYPE *»' — — ' 

DO 305 JX=1fIX 
WRITE ( 5 f 40) JX 

40 FORMAT (' BOUNDARY FLAW SIZE * ' r 13 t ' -(ENTER SIZE IN MICRONS) ') 
READ (5*55) LX 
55 FORMAT (13) 

LR(JX)“LX 
305 CONTINUE 

TYPE * * ' ' 

TYPE Kir' ENTER THE AMOUNT OF OVERLAP (X) DESIRED' 

READ (3*20) N 
20 FORMAT (13) 

TYPE *t ' ' 

Z«0 

CALCULATE OVERLAP FLAW SIZE BOUNDARIES BASED ON 
• OF ORIGINAL BAOUNDARY FLAW SIZES ENTERED* THE 
ORIGINAL SIZES THEMSELVES* AND THE AMOUNT OF OVERLAP 

IEO«I IFIX ( (FLOATI <LR(2) >-FLOATI (LR( 1 ) > )*FLOATI ( 100-N)/100. ) 

IF < IEO.NE.O) GOTO 838 
lEO^l 

838 DO 100 Y=1*IX~1 

DO 200 X=LR(Y)*LR<Y+1)fIE0 

IF (X.EG.LR(l)) GOTO 19? ! DON'T UKIP FIRST BOUNDARY 

IF (X.EQ.LR(Y)) GOTO 200 (DON'T REPEAT BOUNDARIES 

199 Z-Z+l I KEEP TRACK OF HOW MANY BOUNDARY FLAW SIZES 
JR ( Z ) »X (STORE CALCULATED BOUNDARIES IN ARRAY 

200 CONTINUE 
100 CONTINUE 

DO 300 M0=1 t Z (WRITE OVERLAP BOUNDARIES TO TERMINAL 
WRITE ( 5 f 301 ) MOfJR(MO) 

301 FORMAT (' FLAW SIZE * ' f 13 f ' f IX f 13) 

300 CONTINUE 
TYPE #f ' 


THIS ROUTINE PLACES FLAWS IN PROPER INTERVALS 
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300 IQ“0 ! INITIALIZE * OF FLAWS/INTERVAL COUNTER 

KQ“0 I INITIALIZE * OF FLAWS FOUND/ INTERVAL COUNTER 
DO 350 JQ=1,Z IFOR EACH CALCULATED INTERVAL 

IF <JR<JQ>.GT.<.!R<Z)-<LR<2>-LR<i>>)> 00T0 550 
DO 590 NB»1 » IB IFOR EACH FLAW 

IF (LB(NB) .GE. JR< JQ) , AND.LB(NB) • LE , JR( JQ ) + <LR< 2 >-LR( 1 > ) ) 

1 GOTO 505 
GOTO 590 

503 10=10+1 (KEEP TRACK OF • OF FLAWS/ INTERVAL 

IF (KB (NB ) . EQ . 2 ) GOTO 590 I FLAW NOT SEEN ON SLAM 
KQ“K0+1 I KEEP TRACK OF # OF FLAWS SEEN ON SLAM/INTERVAL 

390 CONTINUE 

WRITE (5»591 ) KQ’ IQ. JR( JQ) r JR( JQ) + (LR<2)-LR< 1 ) ) 

591 FORMAT ( IX ’ 13 ’ IX ’ 'OUT OF' . 13’ IX’ 'FLAWS BETWEEN' » 1X» 14’ IX’ 'AND' r 
1I4’1X’ 'MICRONS WERE FOUND ON SLAM') 

TYPE *»' 


AR1=FL0ATI < 10) I CHANGE IQ TO REAL 
AR2 n FL0ATI ( KQ) I CHANGE KQ TO REAL 
AR3“FLOATI(IAR3)/100. 


CALL DIN(AR1 ’AR2rAR3,AR10) { SUBROUTINE FOR PROB. CALC. 


592 WRITE (5 ’593) IAR3.AR10 

593 FORMAT <1X> 'AT'»1X»I2’ ' ’ IX’ 'CONFIDENCE. THE LOWER-BOUND 
1 PROBABILITY OF DETECTION*' . 1X.F4, 2) 

LMX=LMX+1 I KEEP TRACK OF • OF POINTS TO BE PLOTTED 
TYPE #. ' ' 

IF (L32.EQ.0) GOTO 594 (RADIO. PTS. 

XYDATA(2*JQ-1>=50+<JR(JQ>+(LR(2)-LR(1>>>*2 I FLAW SIZE PIXEL COOR. 
XYDAT< 2# JQ-1 > =50+ ( JR < JQ) + ( LR( 2) -LR( 1 ) > ) *2 (THICKEN PLOT LINE 
GOTO 595 I SKIP RADIO. PTS, 

594 XYDATA (2# JQ-1 )»50+NINT < < ( JQ ) + (LR<2)-LR( 1 > ) )*.8) IRADIO. 

XYDAT <2#JQ~1 )=50+NINT < ( JRi JQ)+(LR(2)-LR( 1 ) ) )#,8) IRADIO, 

595 XYDATA(2*JQ)=50+IIFIX(AR10*100.*4. > I P.O.D. PIXEL COOR. 

XYDAT(2*Ja>=49+IIFIX(AR10#100.*4. > (THICKEN PLOT LINE 

PLOT POINTS USING ABOVE COORDINATES (VISIBILITY AID) 

CALL GRFAR ( 1 ’255r0’0»XYDATA(2#JQ-l >-2’XYDATA(2*JQ)-2’4r4> 

PLOT FLAW SIZE RANGE BARS 


IF (L32.EQ.0) GOTO 596 IRADIO. PTS. 

CALL GRFVC ( 1 . 150. 0 . 0 . XYDATA < 2*JQ-1 )-<LR( 2) -LR< 1 ) )$ 2 ? 

1XYDATA < 2* JQ ) . XYDATA < 2* JO-1 ) , XYDATA < 2* JQ ) ) 

CALL GRFVC ( 1 . 150. 0, 0 . XYDATA(2*JQ-1 >-<LR(2)-LR< 1 ) >*2. 

, 1XYDATA(2*JQ)-2’XYDATA(2*JQ“1)-(LR(2)-LR(1) >*2f 
1 XYDATA <2* JO) +2) 

GOTO 597 I SKIP RADIO. PTS. 

596 CALL GRFVC <1 . ISO. 0. Of XYDATA<2*JQ-1 >-NINT< <LR(2)~LR( 1 > )#. 8> » 
1XYDATA(2#JG) r XYDATA ( 2* JQ-1 > ,XYDATA(2*JQ> ) 
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c 

CALL GRFVC ( 1 1 150»0.0.XYDATA<2#JQ-1 )-NINT< <LR<2>-LR( 1 ) )*.8> . 
lXYDATA<2#ja>-2rXYDATA<2*J0~l )~NINT< <LR<2>-LR< 1) >*.0) . 

1 XYDATA ( 2# JQ > +2 ) 

C 

597 10=0 (RESET FLAW COUNTER FOR NEXT INTERVAL 

K0=0 (RESET FLAWS FOUND COUNTER FOR NEXT INTERVAL 
550 CONTINUE 

CALL GRFVL < 1 .255.0.0. XYDATA.LMX) IPLOT OF P.O.D. VS* FLAW SIZE 

CALL ORFVL < 1 .255.0.0.XYDAT.LMX) ITHICKEN PLOT LINE W/XYDAT PTS. 

CALL GRSBFD 

LMX“0 

N99=1000 

TYPE *.' enter: ' 

TYPE 0 - TO GROUP DATA USING ANOTHER METHOD' 

TYPE #»' 1 - TO USE DIFFERENT CONFIDENCE LEVEL' 

TYPE *.' 2 - TO USE DIFFERENT INTERVALS' 

TYPE *»' 3 - TO PLOT • OF FLAWS VS. FL.SZ./TH. SENS.' 

TYFE 4 - TO END PROGRAM' 

READ (5.542) MMM 
542 FORMAT (II) 

TYPE *. ' ' 

IF (MMM.EQ.O) GOTO 611 
IF (MMM.EO.l) GOTO 629 
IF (MMM.EQ.2) GOTO 2178 
IF (MMM. EG. 3) GOTO 2178 
8020 CALL GRSPFD 
CALL GRSEND 
END 

> 
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C THE FOLLOWING SUBROUTINE CALCULATES THE PROBABILITY OF 

C DETECTION OF FLAWS ON THE SLAM (OR OTHER NDE INSTRUMENT I.CMUO'), 

C GIVEN THE SIZE RANGE t CONFIDENCE LEVEL ( AR3) * TOTAL • OF 

C FLAWS ( AR1 ) t AND THE * OF FLAWS FOUND (AR2). 

C 

SUBROUTINE BIN(AR1 t AR2» AR3» ARIO) 

1 IF < AR2 ) 2 » 2 » 4 

2 ARIOSO . 0 

3 RETURN 

4 IF < AR2-AR1 ) 7 » 5 r 5 

5 ARJ0=<1.0-AR3>#*<1.0/AR1) 

6 RETURN 

7 ATT»2.0*AR2 

8 IF (ATT-AR1 )9>9» 12 

9 AR4=AR2-1 . 0 
10 AR5=--i.O 

n GO TO 15 

12 AR4 K AR1- AR2 

13 AR3= 1 . 0-AR3 
L4 AR5*=~1.0 

15 ARIOSO. 5 

16 AHA- 1.0 

17 ARB-0 . 0 

18 AR9*1,0 

19 AR1 1=AR1 

20 AR7=><AR10#*AR8)*< (1 . O-Afil 0 > #* ( AR1-AR8 > > 

21 IF ( AR8-AR4 )22»27»22 

22 AR8=AR8+ 1 • 0 

23 AR9~AR9*AR1 1/ARB 

97 IF (AR9.GT. 10. **30. ) GOTO 98 

24 AR1 1-AR1 1-1 * 0 
96 GOTO 25 

98 AR9=10.**30. 

99 GOTO 24 

25 AR7-AR7+AR9*(AR10**AR8)*< (1 ,0-AR10>**<ARl-AR8> ) 

26 GO TO 21 

27 IF ( AR3-AR7 ) 28 r 28 r 30 

28 AR20=AR10~AR5/<2,0**(AR6+1 .0) > 

29 GO TO 31 

30 AR20=AR10TAR5/(2.0**(AR6+1 .0) ) 

31 CCC=ABS ( AR7-AR3 ) 

32 IF<CCC-0.0001)36f36t33 

33 AR6=AR6+ 1 , 0 

34 AR10=AR20 

35 GO TO 17 

36 lF(ATT-ARl>6i6f37 

37 ARIOSI . 0-AR10 

38 RETURN 

39 END 
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APPENDIX C 


SEEDED SUkFACE VOID INSPECTION DATA 

Figures 13 to 20, generated using the Fortran program described In 
appendix B, show the raw data obtained from the SIAM Inspections of the seeded 
test specimens. Each figure shows the void size distribution plot, that Is, 
the number and size of voids detected and examined for each particular mate- 
rial, specimen thickness, and surface condition. The curves of probability of 
detection (at a 0.95 confidence level) as a function of void diameter (POD) In 
figures 8 to 10 were derived from these data. 
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(a| Schematic diagram o( scanning laser acoustic microscope. 



Ibl Acoustic microscopy technique for ceramic bars. 
Figure - Scanning laser acoustic microscope. 
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Figure 2.- Center and edge optical micrographs 
ot polished, unetched ceramic test specimen 
cross sections showing porosity within 
material. 
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Figure 4, - Optical micrographs ol seeded surface void in sintered silicon 
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change in diameter ol void after polishing. 
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seeded surface void and side-view diagrams 
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Fiqur 1 6. - R«pr fsentitivt surfact profits and corresponding scanning laser acoustic micrographs for s ntered silicon nitride 
spacimtns. 




OWC* . 

OF POC i 



1*1 Diamond ground surl«c* prolil* taken perpendicu (ti Acoustic microqraph ol diamond qround specimen. 

Ur to grinding maiks. 


t igure 6. Concluded. 
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figure 7. Representative surface profiles and '■or espondinq scanninq laser acoustic microqraphs lor sintered silicon carbide 
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Figure 9. - Probability of detection lor voids 
In as-tired sintered silicon carbide spec- 
imens. Probability ol detection calculat- 
ed at 01 95 confidence level. 
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Figure 12. ■ Schematic of steps and plot of re- 
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nominal thickness (39 voids detected/© voids 
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nominal thickness 129 voids detected/40 voids 
total). 
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